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1.0 Introduction 

The ability to accurately predict the thermomechanical deformation response of 
advanced composite materials continues to play an important role in the develop- 
ment of these strategic materials. Analytical models that predict the effective 
behavior of composites are used not only by engineers performing structural anal- 
ysis of large-scale composite components but also by material scientists in devel- 
oping new material systems. For an analytical model to fulfill these two distinct 
functions it must be based on a micromechanics approach which utilizes physi- 
cally based deformation and life constitutive models and allows one to generate 
the average (macro) response of a composite material given the properties of the 
individual constituents and their geometric arrangement. Only then can such a 
model be used by a material scientist to investigate the effect of different deforma- 
tion mechanisms on the overall response of the composite in order to identify the 
appropriate constituents for a given application. However, if a micromechanical 
model is to be used in a large-scale structural analysis it must be 1 ) computation- 
ally efficient, 2) able to generate accurate displacement and stress fields at both 
the macro and the micro level and 3) be compatible with the finite element 
method. Additionally, new advancements in processing and fabrication techniques 
now make it possible to engineer the architectures of these advanced composite 
systems. Full utilization of these emerging manufacturing capablities require the 
development of a computationally efficient micromechanics analysis tool capable 
of accurately predicting the effect of microstructural details on the internal and 
macroscopic behavior of composites. The above mentioned computational effi- 
ciency is required since 1 ) the large number of parameters that must be varied in 
the course of engineering (or designing) composite materials, and 2) the optimiza- 
tion of a material’s microstructure will require the integration of the micromechan- 
ics model with optimization algorithms. From this perspective, analytical 
approaches that produce closed form expressions which describe the effect of a 
material’s internal architecture on the overall material behavior are preferable to 
numerical methods such as the finite element or finite difference schemes. 

A number of models presently exist that can fulfill some aspect of the aforemen- 
tioned tasks. However, there are very few working models that are both computa- 
tionally efficient and sufficiently accurate at the micro- as well as the macro-level. 
One such micromechanics model with the potential of fulfilling bothtasks is the 
method of cells [1] and its generalization [2]. The comprehensive capabilities and 
efficiency of this method has been documented in references [4] and [5]. Conse- 
quently, the recently developed, computationally efficient and comprehensive 
micromechanics analysis code, MAC, who’s predictive capability rests entirely 
upon the fully analytical micromechanics model, herein referred to as the general- 
ized method of cells, GMC, [2 and 3] will now be described. MAC is a versatile 
form of research software that “drives’’ the double or triply periodic micromechan- 
ics constitutive models based upon GMC. GMC is capable of predicting the 
response of both continuous and discontinuous multi-phased composites with an 
arbitrary internal microstructure and reinforcement shape. GMC is a continuum 
based micromechanics model that provides closed-form expressions for the mac- 
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roscopic composite response in terms of the properties, size, shape, distribution, 
and response of the individual constituents or phases that make up the material. 
GMC also utilizes physically based viscoplastic deformation and life models for 
each constituent. Furthermore, expressions relating the internal stress and strain 
fields in the individual constituents in terms of the macroscopically applied 
stresses and strains are also available through strain or stress concentration fac- 
tors. These expressions make possible the investigation of failure processes at 
the microscopic level at each step of an applied load history. Similarily, GMC pro- 
vides the capability of studying the influence of bond strength at the fiber/matrix 
interface which recently has been shown to be an important damage mechanism. 

MAC enhances the basic capabilites of GMC by providing a modular framework 
wherein 1) various thermal, mechanical (stress or strain control) and thermome- 
chanical load histories can be imposed, 2) different integration algorithms may be 
selected, 3) a variety of constituent constitutive models may be utilized and/or 
implemented and 4) a variety of fiber architectures may be easily accessed 
through their corresponding representative volume elements. Figure 1 illustrates 
the basic flow diagram for this modular framework. The capabilities for this version 
of MAC are discussed in section 2, whereas theoretical and background informa- 
tion on the basic capabilities itemized above is given in section 3. Section 4 
describes how one might use MAC and will be the most referred to section in the 
entire manual. Finally, section 5 gives some insights into the future modifications 
planned for MAC. 


3 



MAC M icromechanics A nalysis Code 



Figure t : MAC 
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2.0 Current Capabilities 

In this section the current features/capabilities of MAC are itemized. 

• Load Types: 

- Thermal 

- Mechanical 

- cyclic (stress or strain control) 

- creep (stress control) 

- relaxation (strain control) 

- Thermomechanical 

• Aboudi GMC Models: 

- Double Periodicity Model for continuous reinforcement 

- Triple Periodicity Model for discontinuous reinforcement 

• Run Options: 

- Interactive 

- Batch 

• Graphical Output 

- Up to 5 x-y data plot files may be generated for both macro and micro (sub- 
cell) quantities 

- PATRAN for subcell geometry and color results evaluation, e.g. stress, 
strains, inelastic strains, J 2 , etc. 

• Integration Options: 

- Forward Euler 

- Predictor/Corrector 

• Constitutive Models 

- elastic orthotropic 

- inelastic viscoplastic models: 

- Bodner-Partom 

- Robinson 

- GVIPS 
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3.0 Background 

3.1 Micromechanics Models 

As stated in the introduction, MACs predictive capabilities rest entirely upon the 
fully analytical micromechanics model known as GMC which is capable of predict- 
ing the inelastic response of both continuous (double periodicity) and discontinu- 
ous (triple periodicity) multi-phased composites with an arbitrary internal 
microstructure and reinforcement shape. Prior to describing the available archi- 
tectures (Representative Volume Elements, RVE’s) within MAC as discussed in 
section 4, a brief overview of the theoretical foundation behind the generalized 
method of cells follows. A more complete discussion of the theoretical formulation 
is given in [1] - [3]. 

In the original formulation of the method of cells, a continuously (or discontinu- 
ously) reinforced, unidirectional fibrous composite is modeled as a rectangular, 
double-periodic (or triply-periodic) array of fibers embedded in a matrix phase. 
The periodic character of the assemblage allows one to identify a repeating unit 
cell that can be used as a building block to construct the entire composite. The 
properties of the repeating cell are thus representative of the properties of the 
entire assemblage. The unit cell consists of a single fiber subcell surrounded by 
three matrix subcells for continuous and seven for discontinuous composites, 
hence the name method of cells. The rectangular geometry of the repeating unit 
cell allows one to obtain an approximate solution for the stresses and strains in 
the individual subcells given some macroscopically homogeneous state of strain 
or stress applied to the composite. The approximate solution to the posed bound- 
ary value problem is, in turn, used to determine macroscopic (average) or effec- 
tive properties of the composite and the effective stress-strain response in the 
inelastic region. 

In the generalized method of cells for continuous (or discontinuous) fibrous 
composites, the repeating unit cell can consist of an arbitrary number of phases. 
Hence the generalized method of cells is capable of modeling a multiphase com- 
posite. This generalization extends the modeling capability of the original method 
of cells to include the following: 1 ) inelastic thermomechanical response of mul- 
tiphased metal matrix composite, 2) modeling of various fiber architectures 
(including both shape and packing arrangements), 3) modeling of porosities and 
damage, and 4) the modeling of interfacial regions around inclusions including 
interfacial degradation. 

The basic homogenization approach taken in the micromechanical analysis 
consists essentially of four steps. First, the repeating volume element, RVE, of the 
periodic composite is identified. Second, the macroscopic or average stress and 
strain state in terms of the individual microscopic (subcell) stress and strain states 
is defined. Third, the continuity of tractions and displacements are imposed at the 
boundaries between the constituents. These three steps, in conjunction with 
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micro-equilibrium, establish the relationship between micro (subcell) total, thermal 
and inelastic strains and macro (composite) strains via the relevant concentration 
tensors. In the fourth and final step, the overall macro constitutive equations of the 
composite are determined. These four steps form the basis of the micro-to-macro- 
mechanics analysis which describe the behavior of heterogeneous media. The 
resulting micromechanical analysis establishes the overall (macro) behavior of the 
multi-phase composite and is expressed as a constitutive relation between the 
average stress, strain, thermal, and inelastic strains, in conjunction with the effec- 
tive elastic stiffness tensor. 

That is, 


g = £* (e - e ; - e t ) <EQ i) 

where for the most general case of discontinuous reinforcement with N a by by 
N y number of subcells, the effective elastic stiffness tensor, B , of the composite 
is given by, 


N„ N. N. 


B 


"a '’p 

I 


I W r c 


(<*Py)^ (<*Py) 


(EQ2) 


a = ip = ly = 1 


the composite inelastic strain tensor is defined as, 

e' = £ I (P <oM s',-f' <oW ) 3) 

a = ip = lY=l 

the average thermal strain tensor as, 


f = £ X X d *¥y£ (a?y) (P (aP7) §I-§ r(af5Y) ) (EQ4) 

a = lp = ly = 1 

and e is the uniform applied macro (composite) strain. For the case of continuous 
reinforcements with Np by N y number of subcells, eq. (2) - (4) reduce to the fol- 
lowing: 

5* = I X X (EQ5) 

"'p=lY=l 



-B 


*-i * T 


hi 


X X Vy~ C 


p = ly= 1 


-B 




hi 


X X Vy~ c 


(Pt) 


(Py) 


P=1Y=1 


(P (p7) e^-e /(pY) ) 
(P< PY ) e[-e r( PY ) ) 


(HQ 6) 
(EQ 7) 
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In the above equations matrix notation is employed; where, for example, the 
average stress, g, average applied strain, e, and inelastic subcell strain, e^, vec- 
tors represent, 


O — G 22 ’ ^33’ ^12’ °23’ °13^ 

(HQ 8) 

l = t®ll» ®22* ®33* ®12’ ^23* ^13^ 

(EQ9) 


(EQ 10) 


where the six components of the vector e /(aPY) are arranged as in eq. (9). Similar 
definitions for e Y , e r( “ PY) also exist. Note that the key ingredient in the construc- 
tion of this macro constitutive law is the derivation of the appropriate concentration 
matrices, A (<xPy) and Z) (apY) having the dimensions 6 by 6 and 6 by N a N^N„ 
respectively, at the micro (subcell) level. The definitions of A and D, althouan not 
given here, may be found in references [2] and [5]. Finally, the matrix C (otpY} rep- 
resents the elastic stiffness tensor of each subcell (aPy) and d a , h^, l y the 
respective subcell dimensions (see Fig. 2) wherein, 

N* N y 

<*- S “a *= ;= 

a = 1 3 = 1 y= 1 

Similarly, given the concentration matrices A (apY) and £> (<xPy) , expressions for 
the average strain in each subcell can be constructed, i. e., 


e (aPr) 


= A (afjY) e + D (apY) 



as well as average stress, 


c (otpY) = £<°PTr) [^ a W|+p(“W (e / + e T ) - (e /(aPY) + e r(aPY) )] 


The analytic constitutive law, see eq. 1 , may be readily applied to investigate 
the behavior of various types of composites, given knowledge of the behavior of 
the individual phases. Numerous advantages can be stated regarding the current 
macro/micro constitutive laws as compared to the other numerical micromechani- 
cal approaches in the literature, e.g. the finite element unit cell approach. One 
advantage is that any type of simple or combined loading (multiaxial state of 
stress) can be applied irrespective of whether symmetry exists or not, as well as 
without resorting to different boundary condition application strategies as in the 
case of the finite element unit cell procedure. Another, advantage concerns the 


8 



MAC. M icromechanics A nalysis C ode 


Double 

Periodicity 



/ 



h k 


S 


P=N p 


P = 2 
P-1 




1 

1 


1 



1 1 














'N f 


h 2 

hi 


Y- 1 Y=2 Y=Ny 


A 


h 



Triple 

Periodicity 


a=3 d 3 


a=2 


a=l di 



hi h2 h 3 h4 

P=1 P=2 P=3 P=4 


Figure 2: Subcell Dimension Nomenclature 
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availability of an analytical expression representing the macro elastic-thermo- 
inelastic constitutive law thus ensuring a reduction in memory requirements when 
implementing this formulation into a structural finite element analysis code. Fur- 
thermore, this formulation has been shown to predict accurate macro behavior 
given only a few subcells, within the repeating cell (see references [2], and [4]). 
Whereas, if one employs the finite element unit cell procedure, a significant num- 
ber of finite elements are required within a given repeating unit cell to obtain the 
same level of accuracy as with the present formulation. Consequently, it is possi- 
ble to utilize this formulation to efficiently analyze metal matrix composite struc- 
tures subjected to complex thermomechanical load histories. This is particularly 
important when analyzing realistic structural components, since different loading 
conditions exist throughout the structure, thus necessitating the application of the 
macromechanica! equations repeatedly at these locations. 

3.2 Integration Algorithms 

There are two integration algorithms currently available within MAC. The first is 
the standard, explicit Forward Euler algorithm, which can be expressed as, 

w i+i - w i + A//‘(r | .,w i ) 


where 


w i=y(*i> 
W i+1 =y(r i+1 ) 


and f{t if w,) is the rate of change with respect to time of the vector, y (r 4 ) , i.e. 


/(',•> w f ) = y = 


dyitj 

dt 


The second is a predictor/corrector algorithm which uses: 
1 ) a 4 th order Runge-Kutta starter: 


w i = w i-i + 


(tfj +2K 2 + 2K 3 + K 4 ) 


where 
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K l =Atf(t i _ v w i _ l ) 

A t K, 

K 2 = A + -j,w i _ l + -y ) 
A t 

K 3 = A j + y, VV,--! + y ) 

AT 4 = Atf(t i _ l + At,w i _ 1 +K 3 ) 


2) with an Adams Bashforth four-step predictor: 


w 0 = ctj 

At 


w i = cc 2 


w 2 = a 3 


w 3 = a 4 


wf +1 = w i + 24 [55/(r,-, Wi) -59/(r,._ 1 , w,._ j) +37 /(r,_ 2 > w-_ 2 ) 


-9/(r i _ 3 ,w < ._ 3 ) ] 


and 3) an Adams Moulton four step corrector: 


Ar 


W; 


ui = w i + 24 i’ <+ 1 } + 19/( '<’ ” 5/(r ‘- 1’ w i - 1 } 




where the a’s come from the 4 th order Runge-Kutta starter. Further details may 
be found in [6]. 

It has been found, based on experience, that for relatively rapid monotonic or 
cyclic loadings it may be more efficient to use the Forward Euler integrator since 
the predictor/corrector requires 5 evaluations per step, as shown above. However, 
in the case of creep, relaxation or slow monotonic or cyclic loading histories, sig- 
nificant increases in solution speeds can be obtained using the predictor/corrector 
algorithm with a self-adaptive time step. 
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Finally, within MAC the vector y as used above contains the following macro 
quantities: 


position 

contents 

1 -6 

Macro Total Strain 

7-12 

Macro Stress 

13-18 

Macro Inelastic Strain 

19-30 

(currently empty 
space for 2 6x1 vectors) 

31 -36 

Macro Thermal Strain 

37 

Current Temperature 


and the quantities associated with each subcell are stored sequentially in y, such 
that 


position 

contents 

38-43 

Micro Total Strain 

44-49 

Micro Stress 

50-55 

Micro Inelastic Strain 

56-67 

Micro Internal State Variables 
(space for 2 6x1 vectors) 

68-73 

Micro Thermal Strain 


The above 36 positions are repeated for the total number of (N) subcells thus 
bringing the total length of the y vector to 37+36N. It follows that a second vector 
of similar length contains the corresponding macro and micro rates y. 

3.3 Available Constituent Constitutive Models 

Currently MAC provides four inelastic constitutive models. These models have 
been selected purely based upon the availability of material parameters for the 
materials of interest. However, MAC is designed in a modular fashion thus allow- 
ing the implementation of additional inelastic models. Two of the four available 
models are capable of representing transversely isotropic material behavior, thus 
allowing one to investigate the reinforcement of an anisotropic matrix or the ideal- 
ization of an anisotropic fiber. In all cases a purely elastic response is possible by 
modifying a single material parameter for each model as noted below. 
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3.3.1 Bodner-Partom Model 

Reference: Mechanics of Composite Materials, Jacob Aboudi, Elsevier, 1991 

This model represents the Bodner-Partom viscoplastic model with isotropic 
hardening, Z, and can be used for an initially isotropic metallic material. 

The flow law is given as: 


where 


A = 




A 2 


= -Z 2 

3 Z «// 


(Oil) 


n 


l 

n 


J 2 - 2 S ij S ‘j 
^ij ~ °ij ~ 3 °kk^ij 


The evolution law for isotropic hardening is given as: 

W PL 

Z = m{Z x -Z eff )-^~ 

where Z 0 , Z x and m are inelastic constants and the plastic work rate, W />L , is 
given by; 

* PL = vi 

t 3 t 

Z eff = Z 0 + q jz (r) dx + ( 1 - q) £ r ij\ t r ij (t) 

0 *./ = 1 0 


rijU) = OijO)/ [o w (0o w (r)] 


1/2 


An elastic only response may be obtained by setting the material parameter D 0 to 
zero. 
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3.3.2 Robinson Creep Model 


Reference: NASATM 103172, 1990 


This model represents a transversly isotropic material, wherein the vector of 
direction cosines, d t defines the strong material direction along which no inelastic- 
ity occurs. 

Flow Law: 


e' 

y 


2 F n 


r. 

y 



Evolution Law: 


CL: 


— 1 tn RC m ~ P ^y 

G p 2e 0 o 0 jG 


where 


r , rVl ,( “rV 


in which 


F= 4<t> 

<t, = y 2 -|/ 2 

G = 

°o 



and 


y 2 = j Vi» 


/ = D Z- 
y 7» 


D = dd 
y * y 


^y ^y a v 


/ 2 = X V 


/ = Da 

<7 7* 


s « = »rM 


Special cases involving an isotropic material can be obtained by defining 
= (1/3, 1/3, 1/3) and/or elastic only response by setting o 0 to an extremely 
large number. 
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3.3.3 Robinson Viscoplastic Model 
Reference: NASATM 105819, 1992 

This model represents a transversly isotropic material wherein the vector of 
direction cosines d i defines the perferred material direction. In this model the 
strength of anisotropy is specified by the parameters co and t| ; where co is the 
ratio of the normal longitudinal and transverse yield stress and r\ is the ratio of 
longitudinal and transverse shear strengths. 


Flow Law: 


Evolution Law 



(n r 

2\i y 


x _ H ./ nr m ~ P 

“a - -zi £ ij- RG 


n. 


where 


and 


T ij « V * (D ^ + D A " 2/ ° D y } ' (3D y - 5 y ) 
n ij = a u ~ $ ki a jk + D jk a ki ~ 2 io°ij ) - jtfo (3 D i} - 5 i7 ) 


K 2 


7 1 + Z2 7 2 + 




4 (4co 2 - 1) 


7 3 ]-i 


i r* 


1 ^ 


C = — X | /i + — ^2 + 


2 4(4co 2 — 1) 

G = {G- Gq)Hv [SjjKij] + Gq 


7 l=- 7 2- 7 -Z 7 3 7 2 = 7 — 7 3 7 3 = 7 0 


A = 


Ii.e.. 
2 y V* 


/ = D E 
y y» 


= Wj 


hi 


^y a ‘j 


4 = 


T| 2 - 1 
Tl 2 


C = 


4((0 2 - 1) 
4o) 2 - 1 
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The invariants / lf I 2 , h are the same as those given above but with l.j replaced 
by a t j. Special cases involving an isotropic material and/or elastic only response 
can be obtained by defining co = rj = 1 and/or by setting K r to an extremely large 
number. 
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3.3.4 Generalized Viscoplastic Potential Structure (GV1PS) Model 
Reference: NASATM 106609, 1994 

This model is a fully associative, multiaxial, isothermal, nonlinear kinematic 
hardening viscoplastic model for use with initially isotropic metallic materials. A 
unique aspect of this model is the inclusion of non-linear hardening through the 
use of a compliance operator Q - kl in the evolution law for the back stress. This 
non-linear tensorial operator is significant in that it allows both the flow and evo- 
lutionary laws to be fully associative (and therfore easily integrated) and greatly 
influences the multiaxial response under non-proportional loading paths. 


Flow Law: 


e'. 

y 


^ e. 



if 


F> 0 


where 



2 ./ ./ 
-£• £ 

3 v y 


pF" 

K 


Internal constitutive rate equation 

® ij = L> ijrArs 


Evolution Law: 


if 




-j ft a 3/? a B 0 G^ 

Ar, = <4“ 2jlM i? a ” 

Qn**-, (4-ll«4ll^ m - ) 


o 


if afo „<0 


where 

F = (—-Y) 

K 


y=<i-p7G) 
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Lijrs ~ fCj/rj] 1 


K 2 

3B 0 (l+B lP GP-') 



3B 1 (p-l)G p ~ 2 ~ >| 

k 2 ( 1 + “ 1 (6p - 5) ) a " af; ; 


and 

3 

^2 ~ 2. a ‘j a y — $ij~ a ij 

The special case of an elastic only response maybe obtained by setting k to an 
extremely large value. 
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4.0 Running MAC 

MAC provides the user with the option of either running the code in a interactive 
or batch mode. In the interactive mode, the user is prompted for all necessary 
input through a series of questions. This mode is most useful for beginners or 
when solution of a single problem is desired. Execution in the batch or back- 
ground mode requires the user to construct an input file containing the required 
input as discussed below. This mode is extremely helpful when performing para- 
metric studies in which only a single parameter is being varied. Since the input 
requirements for both modes are identical and the interactive input is self explan- 
atory, only the batch mode input requirements will be presented in detail. First, 
however, a description of the various files required and generated by MAC will be 
given followed by the actual commands required to commence execution of MAC. 

4.1 MAC Input and Output Files 

Input P titput 

mac. key batch, out 

PATRAN files (optional, see section 4.3.11) 
plot data files (see section 4.3.1 9) 

4.2 Interactive/Batch Mode Execution 

Prior to executing MAC the mac.key file must exist. The following is an example 
of the minimum size (3 lines) mac.key file required to run MAC in the interactive 
mode. 

Required input: 

first line: ^problem title 

second line: *JRUN interactive=1 batch=2 

third line: 1 

See Appendix A for a sample mac.key file for running in the batch mode. Each 
control block necessary for running in batch mode is described sequentially in the 
next section. 

MAC has been developed for a Sun IPX workstation running Sun OS 4.1.3. 
Upon linking the associated object modules that comprise MAC, execution is 
begun by simply typing MAC (in upper case) and hitting return. 
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4.3 Input Requirements: Batch Mode 

The following data should appear in the file: mac.key. This mac.key file is the 
job control data file that is read when the batch mode option is used. Since this file 
is read sequentially, the data must appear in the order given here. Also, note that 
the key file must exist and the first two command blocks must alway s be included 
in this key file, i.e. the Header Line and Run Mode blocks as it will define the exe- 
cution option for MAC. 

Note in the following, each block of input data will have its own subsection and 
will typically contain the following information: 

1 ) statement of purpose 

2) declaration of input block 

3) example 

4) notes 


«** Note: Each data block contains the line: 

# comments 

This line is required but the actual text, including the * , are currently not 
used by the code. The present purpose of this line is to make the mac.key file 
easier to organize and read. 
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4.3.1 Header Line: 

Purpose: Define the title of this particular job (80A format). 

«■ Note: the Header Line is limited to one line: 

* problem title 

Example: ^Transverse tensile response of 35% SCS6/Ti--6-4 

4.3.2 Run Mode : 

Purpose: To select either the interactive or batch mode of MAC. 

* run mode comment 

jrun 

j run: 

1 - Interactive mode 

2 - Batch mode 

Example: (to select batch mode) 

*JRUN interactive=1 batch=2 
2 


4.3.3 Output Print Level : 

Purpose: To control the output generated 

* output comment 
nplvl 

nplvl: 

-1 = print out macro stiffness matrix, engineering constants and macro ther- 
mal expansion coefficients 
0 = minimal output (most commonly used) 

3 = basic program trace 
5 = program execution trace and all array data 

WARNING: this generates a very large output file 

«* Note: When using the batch mode, all output is written to the file, batch.out. For 
the interactive mode, all output is written to the screen. 

Example: (minimal print out) 

*NPLVL 

0 
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4.3.4 Load Type : 

Purpose: To select load type 

* load type comment 

nsel 

nsel: 

1 = Thermal Load 

2 = Mechanical Load 

3 = Thermomechanical Load 

Example: (run mechanical load) 

*NSEL 

2 

Note: This Is where the data would begin if the new problem option is chosen in 
section 4.3.20 

4.3.5 Load Component : 

«*• Note: This data block is NOT TO BE USED if NSEL = 1 (thermal load) 

Purpose: To select the load component (refer to figure 3 for coordinate system to be 
activated) 

# load component comment 

lop 

lop: 

1 = axial load in 1 -direction 

2 = axial load in 2-direction 

3 = axial load in 3-direction 

4 = shear load 1 2-direction 

5 = shear load 23-direction 

6 = shear load 13-direction 

Example: (axial load in 1 -direction) 

*LOP 

1 
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4.3.6 Load Control: 

«sr Note: This block is only required if NSEL = 2 or 3 

Purpose: Select type of load control for mechanical load 

* load control comment 

iopt2 

iopt2: 

1 = Strain control 

2 = Stress control 

Example: (strain control) 

*IOPT2 

1 
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For Thermal Load (NSEL = 1 ): 

* number of time and temperature points comment 

nptt 

* time points comment 
h> *2> tnptt 

* temperature points comment 

T~1, Tg T n pff 

For Mechanical Load (NSEL = 2): 

* number of time and load points comment 
nptw 

* time points comment 

*1> *2> •••» tnptw 

* load points comment 

Lit 1-2, ...» l-npfw 

For Thermomechanical Load (NSEL = 3): 

* number of time and load points comment 

nptw 

* time points comment 

h> h> •••» tpptw 

* load points comment 

!—f, 1-2 , ...» Lnptw 

* number or time and temperature points comment 

nptt 

* time points comment 
f 2> ■> tnptt 

* temperature points comment 

Tf, T 2 , .... T n ptf 


npt, nptw, nptt: 
t"f, t2, •••« 

L l> L 2 , ... : 

Tf, T2, — 


number of points on load and/or temperature curve(s) 
time values 
load curve values 
temperature curve values 


Note: For the Thermomechanical Load t nptw = t npW and both curves must have t 1 = 0. 
But the number and time value of the data points in-between maybe different, see fig- 
ure 4. 
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Examples: 

For thermal load only analysis (NSEL = 1 ), using 4 time and temperature points: 

* NPT 
4 

*TIME 
0..1.. 2., 3. 

* TEMPERATURE 
0., 100., 200. ,300. 

For mechanical load only analysis (NSEL = 2), using 4 time and load points: 
*NPT 
4 

* TIME 
0.,1.,2.,3. 

* LOAD 

0., 10., 20., 30. 

For thermomechanical load analysis (NSEL = 3), using 4 points for both load and 
temperature data: 

*NPTW 

4 

* TIME 
0.,1.,2.,3. 

* LOAD 

0., 10., 20., 30. 

*NPTT 

4 

*TIME 

0.,1.,2.,3. 

^TEMPERATURE 
0., 100., 200., 300. 


i®- Note: Units must be consistent with selected material model parameter units (4.3.14) 
See figure 4 for sketch of load history. 
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Figure 4: Load History Specification 
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4.3.8 Reference Temperature: 

Purpose: Define initial (starting) temperature 

* reference temperature comment 

tref 

tref: 

-- temperature 

Example: (set reference temperature at 70.) 

*TREF 

70. 

m- Note: only important in thermal or thermomechanical load case when using secant 
values of coefficient of thermal expansion. 

< 3 * Note: Units must be consistent with selected material model parameter units 
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4.3.9 Micromechanics Model Identification: 

Purpose: Select desired GMC micromechanics model 

♦ micromechanics model comment 
modid 

modid: 

1 = double periodicity (continuous reinforcement) 

2 = triple periodicity (discontinuous reinforcement) 

Example: (select double periodicity model) 

♦MODID 

1 

«*■ Note: This is where the data would begin if a new case option is chosen in sec- 
tion 4.3.20 
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4.3.10 PATRAN Output: 

Purpose: To select PATRAN output file generation. 

* PATRAN comment 
patran 
prefix 

patran: 

y = yes, generate PATRAN files 
n = no, do not generate PATRAN files 
prefix : 

this will be the filename prefix that will be assigned to all PATRAN files. 

Note: This Is only required if patran is yes 

Example 1: (select PATRAN output, use file prefix “runl”) 

♦PATRAN 

y 

runl 

Example 2: (do not select PATRAN output) 

♦ PATRAN 
n 

**■ Note: If this option is chosen, MAC generates the following PATRAN files: 

1) prefix.p atgeo: contains “geometry” information of the RVE. 

(this file is written in PATRAN 2.5 neutral file format in which the subcells are 
treated as “elements”) 

2) prefix. patstr: contains stress quantities 

(file is formatted as a PATRAN 2.5 element results file) 

3) prefix. patepsin: contains inelastic strain quantities 

4) prefix. patepsto: contains total strain quantities 

For Example 1, the following files would be produced: 

runl.patgeo 
runl. patstr 
runl. patepsin 
runl .patepsto 
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4.3.11 Integrator Identification 

Purpose: Select type of integration scheme 

¥ integrator comment 

ntf 

* initial time step comment 

istp 

* error tolerance comment 

errtol 

ntf: 

1 = Forward Euler method 

2 = Predictor/Corrector method 

istp: 

= initial time step 

errtol: 

= error tolerance for predictor/corrector 
Note: errtol is only required when using predictor/corrector (ntf = 2) 
Suggested errtol = 0.1 


Example 1: (select predictor/corrector, with initial time step = 0.001 and error toler- 
ance = 0.01) 

*NTF 

2 

*ISTP 

0.001 

* ERRTOL 
0.01 

Example 2: (select forward euler, time step = 0.001) 

*NTF 

1 

*ISTP 

0.001 
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4.3.12 Constituent Material Model Identification: 

Purpose: To select the model for the fiber and matrix constituents 

* material model comment 
ncmd 

ncmd: 

1 = Bodner-Partom Model 

2 = Robinson Creep Model 

3 = Robinson Non-normalized Viscoplastic Model 

4 = Generalized Viscoplastic Potential Structure (GVIPS) Model 

Example: (select Robinson Non-normalized Viscoplastic Model) 

*NCMD 

3 

is- Note: Currently all phases, e.g. both the fiber and matrix, constituents, must use the 
same material model. 

■s* Note: See Section 3.3 for a mathematical description of each material model. 
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4.3.13 Number of different Materials: 

Purpose: Define the total number of different materials 

«* Note: see below for the proper material constants required. 


* number of materials comment 
nmat 

* material # 1 comment 
mdata 

* database option comment 
idb 


* material #nmat comment 
mdata 

* database option comment 
idb 



repeat each block nmat 
times 


nmat - number of materials 

mdata - material identification letter selected from material database below 


MAC Isothermal Materials Database 


Model (4.3.13) 

Material 

Units 

mdata 

Bodner-Partom 

Boron 

MPa, sec 

A 

ncmd= 1 

Aluminum 

11 u 

B 


SCS-6 

U U 

C 


Ti-15-3 

u u 

D 


Ti-6-4 

u u 

E 

Robinson Creep 

Tungsten 

ksi, hr 

A 

ncmd= 2 

Kanthal 

u u 

B 

Robinson Viscoplastic 

Tungsten 

ksi, hr 

A 

ncmd = 3 

Kanthal 

u u 

B 


FeCrAlY 

u ii 

C 


Weak interface 

u u 

D 

GVIPS 

TIMETAL21S 

ksi, sec 

A 

ncmd = 4 

SCS-6 

u u 

B 


Weak interface 

u u 

C 
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idb - read material constants from the database (y) or read from mac.key 
file (n) 

ns* Note: if idb is n (no) then the following four lines must also appear: 

* dummy keyword 
elastic constants 

* dummy keyword 
viscoplastic constants 

Example: (select 2 materials, material 1 : Boron, read from database , material 2: Alu- 
minum, enter properites) 

*NMAT 

2 

*MAT1 

A 

*IDB 

y 

*MAT2 

B 

*IDB 

n 

^ELASTIC 

E, v, G, a (see below for actual constants required) 

^INELASTIC 

a, b ,c ,d, ... (see below for actual parameters required) 

Required Format for Isothermal Material Constants: 

Bodner-Partom: ncmd= 1 


Elastic: 

Inelastic: 

E, v 

D 0 , Z 0 , Zi , m, n.q 

Robinson Creep: ncmd= 2 

Elastic: 

Inelastic: 

E-|, E 2i v r v 2 , G 12 , a 2 

Oq, e 0 , n, p, m, R, H 

Robinson Viscoplastic: ncmd 

= 3 

Elastic: 

Inelastic: 

Ei, E 2i Vj, v 2 , G 12 , a,, a 2 

n, m, p, Kj, p, R, H, G 0 , q, to 

GVIPS: ncmd =4 

Elastic: 

Inelastic: 

E, v 

: ji, k, ^? a , /? K , b v P* w, z 0 
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Non-isothermal Material Constants 

The Following materials are available for a non-isothermal analysis: 

MAC Nonisothemal Materials Database 


Model (4.3.13) 

Material 

Units 

mdata 

Bodner-Partom 

Graphite 

MPa, sec, 0 C 

A 

ncmd = 1 

Aluminum 

u u 

B 


User Defined 

u u 

C 


Required Format for Non-isothermal Material Constants: 
ntpts 

T-j, T2, T n tpts 


Ej , , ... Ef 


ntpts 


V y V Op y « e t V ^ 

1 1 1 2 1 ntjHM 


Gy 9 Gy^ ...G n 


ntpts 


1 ntpts 


(Xy y (Xy y • * • Cty 

d o t > d o t ’ -Do, 
2 0r > Z 0r , •••Z 0 

7 1 7 2 'ii/pii 

7 1 7 2 7 «/pu 

m T , m r y ... m T 
7 1 7 2 1 

» • • • 

7 1 7 2 1 ntpts 


<7r> Qts —Qt 

Note: Currently, only the Bodner-Partom Model has nonisothermal constants. 
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4.3.14 RVE Identification : 

Purpose: Select RVE representing desired fiber packing arrangement/architecture 

* RVE comment 

idp 

idp: 


• 1 = 


• 2 = 


• 3 = 


• 4 = 


• 6 = 


Square Fiber, Square Pack (original 4-cell model) 
vy< i/ (i + A) 2 
RVE shown in Fig. 5 

Square Fiber, Triangular (hexagonal) Pack 

V f < 0.86602/(1+ A) 2 

RVE shown in Fig. 5 
Square Fiber, Square Diagonal Pack 

vy<0.5/ (1 + A) 2 

RVE shown in Fig. 5 
Cross Shaped Fiber, Square Pack 

V f < 1 -4 (;ca) 2 

RVE shown in Fig. 5 

Circular Fiber Approximation Rectangular or Square Pack 



R> 1.0 
R< 1.0 
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9 = 


11 = 


• 99 = 


Two Different Size Square Fibers, Rectangular or Square Pack 

1 


v 


R 


Rf ^ 2 

(l+A^ + d+Aj)^ 
fl 


v /, £ 


2 M: 




RVE shown in Fig. 5 
Square Fiber, Rectangular Pack 

1 


V f < = 

f R( 1+A) 2 

V f < ^ 2 

f (1+A) 2 


if tf>1.0 


1.0 


RVE shown in Fig. 5 
User Defined RVE 

Example of RVE representing random packing shown in Fig. 6 
Required input shown in Appendix D. 


Example: (select RVE idp=1 ) 

♦RVE 

1 

Note: R, R f . R f , Vp A = cper, Aj = cperl and A 2 = cper2 are defined in sec. 4.3.16 

ra* Note: The images of the RVE’s shown in Figs.5 and 6 were generated using the PAT- 
RAN option (section 4.3.11) available in MAC. 
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IDP = 11 

Figure 5: RVE’s Available in MAC 
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IDP = 99 


Figure 6: User defined RVE 
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4.3.15 Fiber/Matrix Interface Option: 

♦ fiber/matrix interface comment 

inter 

inter : 

1 - with interface 

2 - no interface 

Example: (select interface option) 
♦INTERFACE 
2 
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4.3.16 RVE Data : 

Purpose: Enter required RVE data. 


Without Interface: (inter=2) 

• For IDP = 1, 2, 3 

* RVE comment 

vf 

• For IDP = 4 

* volume fraction comment 

vf 

* xa distance comment 

xa 

• For IDP = 6 

* volume fraction comment 

vf 

* side ratio comment 

R 

• For IDP = 9 

* volume fraction 1 comment 

vfl 

* fiber radius 1 comment 

radl 

* volume fraction 2 comment 

vf2 

* fiber radius 2 comment 

rad2 

* side ratio comment 

R 

• For IDP = 11 

* volume fraction 1 comment 

vfl 

* fiber radius 1 comment 

radl 

* side ratio comment 

R 

• For IDP = 99 

* RVE dimensions comment 

nb,ng 

* h dimensions comment 

h2 h ng 

* I dimensions comment 

h» b ^nb 
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With Interface: (inter=1 ) 

• For IDP = 1, 2, 3 

* fiber volume ratio comment 

vf 

* fiber radius 1 comment 

radl 

* interface thickness comment 

cper 

• For IDP = 4 

«■ Currently NOT Available 

• For IDP = 6 

is* Currently NOT Available 

• For IDP = 9 

* fiber 1 volume ratio comment 

vfl 

* fiber radius 1 comment 

radl 

* interface 1 thickness comment 

cperl 

* fiber 2 volume ratio comment 

vf2 

* fiber radius 2 comment 

rad2 

* interface 2 thickness comment 

cper2 

* side ratio comment 

R 

• IDP = 11 

* fiber 1 volume ratio comment 

vfl 

fiber radius 1 comment 

radl 

* side ratio comment 

R 

% interface thickness comment 

cper 

• For IDP = 99 

* RVE dimensions comment 

nb.ng 

* h dimensions comment 
h-(, h2 h n g 

* / dimensions comment 
^2> I nb 
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where: 

vf, vfl, vf2 = the fiber volume ratios 
rad, radl, rad2 = fiber radii 

cper, cperl, cper2 = ratios of interface thickness to fiber radius 
xa = length of the cross 

R = X/Y which defines the ratio of distances between fiber within a “ply” and those 
between a “ply” (see figure on following page) 


Hybrid Composite RVE 
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Cross Shaped Fiber Distance xa 



Fiber Spacing Ratio, R=X/Y 
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4.3.17 Plot Point Information : 

Purpose: Specify the frequency at which data will be written to output files for both X- 
Y data and PATRAN if applicable. This provides the user with the ability to control the 
size of the output files since PATRAN files can become large, depending on the prob- 
lem. 

* plot point frequency comment 
npmax 

Example: (print out every fifth data point) 

*NPMAX 

5 
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4.3.18 Curve Data : 

Purpose: Select plot variables for x and y axes and plot data filename 



* curve 1 comment - (curve 1) 

* x variable comment - (x variable) 
maidx 

* y variable comment - (y variable) 
maidy 

* filename comment n- (filename) 
tname 

* continuation comment - (more) 
more? 


maidx and ma/'dy variable options: 


1 - e u 

? ■ °n 

13 * e n 

19- £** 

22 - Total Time 

2 - e 22 

8 -Ojj 

14-42 

20 - £-[2 

23 - Creep Time 

3 - e 33 

9-o 33 

15- 4 3 

21 - 4 

24 - Temperature 

4 -t 12 

10-i 12 

15-42 



5- Y 23 

11 - t,. 

17-43 



6 -t 13 

12 - ’'.I 

18-43 




tname: 

— = name of plot file 

isr Note: The file(s) generated will be of the form tname_ mac.data. If the user desires to 
to use the same tname for files 2-5, a double quote, “, is entered for tname. 
Those files then will have the form: 

tname 2_mac.data 
tname 3_mac.data 
etc. 

(see the example for more details) 

more?: 

y = generate another data file 
n = finished selecting data files 

*3“ Note: Currently a maximum of § curves maybe specified per problem 
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EXAMPLE: 


* CURVE 1 
*X 

1 

*Y 

7 

* FILENAME 

stress file created: stress_mac.data 

* MORE? 

y 

* CURVE 2 
*X 

22 

*Y 

13 

* FILENAME 

“ file created: stress2_mac.data 

* MORE? 

y 

* CURVE 3 
*X 

23 

*Y 

13 

* FILENAME 

plot file created: plot_mac.data 

* MORE? 
n 
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Specifying Curve Data for Micr o (subcelll Quantities: 

* micro option comment 

micro 

micro: 

y = generate subcell data files 
n = no subcell data requested 

if micro is y (yes), then enter the following data lines: 

* number of subcells comment 
nucell 

* curve 1 comment - (curve 1) 

* x variable comment - (x variable) 

miidx 

* y variable comment - (y variable) 

miidy 

* filename comment - (filename) 
tname2 

* dummy keyword 

more? 


miidx and miidy: 


1 - e n 

7 ■ °n 


19- *ii 

25 - W u 

31 - Total Time 

2 ' e 22 

8 -Ojj 

14 -£<2 

20 - <D 22 

26-^ 

32 - Creep Time 

£ 

CO 

1 

CO 

9-o 33 

15-4 

21 - 0 33 

27 - Y* 

33 - Temperature 

4 -?12 

10-t 12 

16-4 

22 - 0 12 

28 - Y 12 


5 - Y 23 

11 - x^ 

17-4 

23-0^ 

29-^ 


6 - Yo 

12 - t 13 

18-4 

24 - 0 13 

30 - ^13 



The quantities <$ and 'F are the possible internal state variables (constitutive model 
dependent). 
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tname2: 

— = name of plot file 

«■ Note: The file(s) generated will be of the form fname2_mic.data. If the user desires to 
to use the same tnam&2 for files 2-5, a double quote, ", is entered for tname 2. 
Those files then will have the form: 

ti 7ame22_mic.data 
tname 23_mic.data 
etc. 

(see the example for more details) 

nucell: 

— = subcell number 
more?: 

y = to enter another subcell number 
n = finished selecting subcells 


EXAMPLE: * MICRO? 

y 

* NUCELL 
1 

* CURVE 1 
*X 

1 

*Y 

7 

* FILENAME 

cell file created: cell_mic.data 

* MORE? 
n 
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4.3.19 Multiple Run Option : 

Purpose: Specify if a new “case” or “problem” is to be run 

* case comment 

ncase 

* problem comment 

nprob 

Example: ENCASE 

n 

*NPROB 

n 


Note: Here, a new “case” means to run the problem with the same load type and his- 
tory but with a new material model , new mate rial constants, and/or RVE. 

A new “problem" is similar to the above except that a new load type and history 
may also be specified. 

This option allows a sequence of runs to be made with different load types, 
load histories, material models, material constants, and packing arrangements 
all contained within a single MAC.key file 

Note: if ncase is y then the data for the next case must follow, for example, 


first problem data 


example of an 
additional case 
desired 

new case data 

start at section 4.3.10 


* NCASE? 
n 

* NPROB? 

n end of mac.key file 


* NCASE? 

y 
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«* Note: if ncase is n but nprob is y then the following format is used, 



* NCASE? 
n 


* NPROB? 

y 



* NCASE? 


n 

* NPROB? 
n 


first problem data 


new problem data 
start at section 4.3.4 


end of mac. key file 


example of an 
additional problem 
desired 


Refer to Appendix B for a detailed mac.key file example. 
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5.0 Conclusion/Future Modifications 

A computationally efficient, user-friendly, comprehensive, micromechanics analysis 
tool, MAC, has been presented that admits physically based viscoplastic deformation 
and life models, can analyze multiphased materials of interest in advanced propulsion 
systems, and can assist both the material scientist and structural analyst in developing, 
designing and analyzing strategic materials. However, the development of this tool is far 
from complete. A number of future enhancements that are planned and currently under- 
way include: 

1) Provide a convenient link for user definable constitutive material models. 

2) Provide an implicit integration algorithm to improve computational efficiency. 

3) The ability for each subcell to utilize a different material model. 

4) The incorporation of damage evolution laws and failure criteria so as to provide 
life estimates. 

5) The incorporation of lamination theory to allow analysis of laminates. 
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6.0 Appendix A 


Sample mac.kev File For A Mecha nical Load Problem 


The following example is used to explain the control blocks in more detail. 
Problem Summary : 


Load Type: 

Load Component: 
Load History: 

Load Control: 

Load History Data: 


Mechanical 

22-direction (transverse to fiber) 

Cyclic 

Strain 

e = 8.333-xlO -4 , 


Micromechanics model: 
Fiber Packing Arrangment: 
Integration Algorithm: 
Constituent Material Model: 
Constituents: 


= 0-015, 

£ min = 0 ‘ 

^initial = 0-0000024 
Double Periodicity 

Hexagonal Pack at 35% fiber volume ratio 

Predictor/Corrector 

GVIPS 

Fiber: SCS-6 
Matrix: TIMETAL21S 

Interface: fictitious weak interface for TIMETAL21 S 
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♦FOR MECHANICAL LOAD RUN 


♦JR UN INTERACTIVE^ BATCH=2 
2 


♦NPLVL 

-1 


♦NSEL 

2 


♦LOP 

2 


♦IOPT5 

1 


♦NPTW 

2 

♦TIMES 

0..18. 

♦LOAD 

0..0.015 


♦TREF 

650. 


♦MODID 

1 


♦NTF 

2 

♦ISTP 

0.0000024 

♦ERRTOL 

0.1E-2 


♦PATRAN 

n 


♦NCMD 

4 


♦NMT 

3 


Header Line 
Batch Control Mode 

Print Level 

Load Type: Mechanical 

Load component in 
22-direction 

Strain Control 

2 Data points 
Time Data 
Load Data 

Reference Temperature 

Double Periodicity 

Predictor/Corrector 
Initial Step = 0.0000024 
Error Tolerance 

No PATRAN Output 

GVIPS Model 

3 materials 
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•MATT 

B 

*IDB 

y 


♦MAT2 

A 

*IDB 

y 


♦MAT3 

C 

*IDB 

n 

♦ELASTIC 

11700..11700., 0.365, 0.365,4287.5,1., 1. 

♦VISCOPLASTIC 

0.8E-8,0.1,0.1E-5,0.,0.85E-3,0.05,1.,1., 1.3.3,1.8,1.35,1., 0.01 


♦IDP 

2 


♦INTERFACE 

1 


♦VF 

0.35 

♦RAD 

0.07 

♦CPER 

0.1 


♦FIBID 

1 

♦INTID 

3 

♦MATID 

2 


♦NPMAX 

5 


♦CURVE 1 
♦X 
2 

♦Y 

8 


Material 1 : SCS-6 
read from database 

Material 2: TIMETAL21S 
read from database 

Marerial 3: Weak Interface 
Read constants from key file 
Elastic constants 
Viscoplastic constants 

RVE: Hexagonal Pack 

Include fiber/matrix 
Interface 

Volume Fraction 
Fiber Radius 
Interface 

Thickness/Fiber Radius 

Fiber Material id. = 1 
Interface Material id. = 3 
Matrix Material id. = 2 

Print requested data every 
5 Increments 

Curve 1 data: 
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♦FILENAME 

SEtranwhex35 


♦MORE? 

n 

♦MICRO? 

n 


♦NCASE? 

n 

♦NPROB? 

n 


Filename for Curve 1 

No more macro quantities 
No micro quantities 

No new cases 
No new problems 


The following figure was obtained from the x-y plot data file produced by the present 
example. 



ef Note: It is recommended that a new user construct a mac.key file using the data given 
in this appendix and then check to see if the same result plot is obtained. 
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7.0 Appendix B 


Sample of A Series of Runs Contained Within One mac.kev File 

♦FOR MECHANICAL LOAD RUN 
♦JRUN INTER ACTIVE=1 BATCH=2 
2 

♦NPLVL 

-1 

♦NSEL 

2 

♦LOP 

2 

♦IOPT5 

1 

♦NPTW 

2 

♦TIMES 

0..12.0 

♦LOAD 

0..0.01 

♦TREF 

873. 

♦MODID 

1 

♦NTF 

2 

♦ISTP 

0.0048 

♦ERRTOL 

0.1E-2 

♦PATRAN 

n 

♦NCMD 

4 

♦NMT 

3 

♦MAT1 

B 

♦IDB 

y 

♦MAT2 

A 

♦IDB 

y 

♦MAT3 

C 

♦IDB 

n 

♦ELASTIC 

1170.,1170.,0.365,0.365,4285.7,1.0,1.0 
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★VISCOPLASTIC 

0.8D-7,0.1,0.1D-5,0.,0.85E-4, 0.05, l.,l.,l.,3.3, 1.8,1.35,1. ,0.01 
*IDP 
1 

♦INTERFACE 

1 

♦VF 

0.35 

♦RADI 

0.07 

♦CPER 

0.1 

♦FIBID 

1 

♦INTID 

3 

♦MATID 

2 

♦NPMAX 

1 

♦CURVE 1 
♦X 
2 

♦Y 

8 

♦FILENAME 

test 

♦MORE? 

n 

♦MICRO? 

n 1 

♦NCASE? 

y _ 

1 

♦NTF 

2 

♦ISTP 

0.0048 

♦ERRTOL 

0.1E-2 

♦PAIR AN 

n 

♦NCMD 

4 

♦NMT 

3 

♦MAT1 

B 

♦IDB 

y 

♦MAT2 

A 


Select new case option 

Beginning of new case data 
Note that it begins at 
section 4.3.10 
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*IDB 

y 

*MAT3 

C 

*IDB 

n 

♦ELASTIC 

11700.,11700., 0.365,0.365,4285.7,10,1.0 
♦VISCOPLASTIC 

0.8D-7,1.,0.1D-5,0.,0.85E-4,0.05,1.,1.,1.,3 .3,1.8, 1.35,1.,0.01 
♦IDP 
1 

♦INTERFACE 

1 

♦VF 

0.35 

♦RAD 

0.07 

♦CPER 

0.1 

♦FIBID 

1 

♦INTID 

3 

♦MATID 

2 

♦NPMAX 

1 

♦CURVE 1 

♦X 

2 

♦Y 

8 

♦FILENAME 

test2 

♦MORE? 

n 

♦MICRO? 

n 

♦NCASE? 

n 

♦NPROB? 

y 

♦NSEL y/////////////////////////////////////////////////////////////////. 
2 

♦LOP 

2 

♦IOPT5 

1 

♦NPTW 

2 

♦TIMES 

0..12.0 


Skip new case option 

Select new problem option 

Beginning of new problem 
data 

Note that it begins at 
section 4.3.4 
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♦LOAD 

0..0.01 

♦TREF 

873. 

♦MODID 

1 

♦NTF 

2 

♦ISTP 

0.0024 

♦ERRTOL 

0.10E-2 

♦PATRAN 

n 

♦NCMD 

4 

*NMT 

3 

♦MAT1 

B 

♦IDB 

y 

♦MAT2 

A 

♦IDB 

y 

♦MAT3 

A 

♦IDB 

n 

♦ELASTIC 

1170.,1170.,0.365,0.365,4285.7,1.0,1.0 

♦VISCOPLASTIC 

0.8D-7,0.1,0.1D-5,0.,0.85E-4,0.05,1.,1.,1.,3.3, 1.8,1.35,1., 0.01 
♦IDP 
1 

♦INTERFACE 

1 

*VF 

0.35 

♦RAD 

0.07 

♦CPER 

0.1 

♦FIBID 

1 

♦INTID 

3 

♦MATID 

2 

♦NPMAX 

1 

♦CURVE 1 
♦X 
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2 

*Y 

8 

•FILENAME 

test3 

•MORE? 

n 

•MICRO? 

n 

End of mac.key file 



62 



MAO M icromechanics 4 naiysis C ode 


8.0 Appendix C 

Sample mac.kev File For A Thermal L oad Problem 
The following example is used to explain the control blocks in more detail. 


Problem Summary : 
Load Type: 

Load History Data: 


Micromechanics model: 
Fiber Packing Arrangment: 
Integration Algorithm: 
Constituent Material Model: 
Constituents: 


Thermal 
T = 0.01, 


T = 371 1 

1 max - ,/ i-1 ’ 


T- =21.1 

min 


= 11505 


Double Periodicity 

Square Pack at 35% fiber volume ratio 

Predictor/Corrector 

Bodner-Partom 


Fiber: Graphite 
Matrix: Aluminum 


Note: This problem is taken from reference 1 ., pg.238 
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•FOR MECHANICAL LOAD RUNT 


♦JRUN INTER ACTIVE= 1 BATCH=2 
2 


♦NPLVL 

0 


♦NSEL 

1 


*NPTT 

3 

♦TIMES 

0.,35010.,70020. 

♦TEMPS 

371.1,21.1,371.1 


♦TREF 

371.1 


♦MODID 

1 


♦NTF 

2 

♦ISTPT 

17.505 

♦ERRTOL 

10E-2 


♦PATRAN 

n 


♦NCMD 

1 


♦NMT 

2 


♦MAT1 

A 

♦IDB 

y 


♦MAT2 

B 

♦IDB 

y 


Number of Data Points 
Time Data Points 
Temperature Points 


Predictor/Corrector 
Initial step size 
Error Tolerance 
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The following figure was obtained from the x-y plot data file produced by the present 
example. 


xlO 1 



«■ Note: It is recommended that a new user construct a mac.key file using the data given 
in this appendix and then check to see if the same result plot is obtained. 
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9.0 Appendix D 

Sample mac.kev File For Rand om Composite 

♦FOR MECHANICAL LOAD RUN 
♦JRUN INTERACTIVE=1 BATCH=2 
2 

♦NPLVL 

-1 

♦NSEL 

2 

♦LOP 

1 

♦IOPT5 

2 

♦NPTW 

2 

♦TIMES 

0..0.04 

♦LOAD 

0..4.0 

♦TREF 

873. 

♦MODID 

1 

♦NTF 

1 

♦ISTP 
0.004 
♦PAIR AN 
n 

♦NCMD 

3 

♦NMT 

2 

♦MAT1 

B 

♦IDB 

y 

♦MAT2 

D 

♦IDB 

y 

♦IDP 

99 

♦INTERFACE 

2 

♦NBNG 

14,14 

*XH 

L, 

1 , 

L, 

L, 
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1., 

1., 

In 

In 

In 

In 

In 

In 

In 

In 

*XL 

♦MATNUM 

laaaaxkAOXi'i&i 

12.UW. 1.1,2A2A2,1 

TziAaazazzuw 

2A2A2A2.1 .1A1.1A2 
2Zl,U^,l,UAW^ 
1.2.1.1A2A2A2.M.2.1 
1^,2ZU,2^,1,U,1 
22,2,2,2,1,1,2,2,2,2,2.2,2 
2,1.1.2A2A2.1,1A1.U 

2.1, U,l,UAl,U,l,U 

2222.1.1.1, UA2A22 
l.UWAl.Uttl.U 

uwwwu; 

2A2A2A2A2AUA2 

♦NPMAX 

1 

♦CURVE 1 
♦X 
1 

*Y 

7 

♦FILENAME 

test 

♦MORE? 

n 

♦MICRO? 

n 

♦NCASE? 

n 

♦NPROB? 

n 
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